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ABSTRACT 



, Received accepted 

^ ■ 

Context. X-ray observations performed with the Rontgen Satellite (ROSAT) led to the discovery of a group (seven to date) 
^ of X-ray dim and radio-silent middle-aged isolated neutron stars (a.k.a. XDINSs), which are characterised by pure blackbody 
^ . spectra (kT « 40 — 100 eV), long X-ray pulsations (P = 3 — 12 s), and appear to be endowed with relatively high magnetic 
CZ3 ' fields, [B ^ 10"-10" G). Optical observations of XDINSs are important, together with the X-ray ones, to study the cooling of 
the neutron star surface and to investigate the relation betwen XDINSs and other isolated neutron star classes. RBS 1774 is 
one of the few XDINSs with a candidate optical counterpart, which we discovered with the Ven/ Large Telescope ( VLT). 
Aims. We aim at constraining the optical spectrum of RBS 1774, for which only two B-band flux measurements are available, 
^ and to determine whether its optical emission is either of thermal or of non-thermal origin. 
^\ Methods. We performed deep observations of RBS 1774 in the R band with the VLT to disentangle a non-thermal power-law 
0^ ' spectrum from a Rayleigh- Jeans, whose contributions are expected to be very much different in the red part of the spectrum. 
, Results. We did not detect the RBS 1774 candidate counterpart down to a 3a limiting magnitude of i? ~ 27. The constraint 
on its colour, (B — R) ^ 0.6, rules out that it is a background object, positionally coincident with the X-ray source. Our 
fT^ , R-band upper limit is consistent with the extrapolation of the B-band flux (assuming a 3a uncertainty) for a set of power-laws 
. Fi, oc v'" with spectral indeces a < 0.07. If the optical spectrum of RBS 1774 were non-thermal, its power-law slope would be 
very much unlike those of all isolated neutron stars with non-thermal optical emission, suggesting that it is most likely thermal. 
For instance, a Rayleigh- Jeans with temperature To = 11 eV, for an optically emitting radius ro ~ 15 km and a source distance 
d = 150 pc, would be consistent with the optical measurements. The implied low distance is compatible with the 0.04 X-ray 
-1 — I pulsed fraction if either the star spin axis is nearly aligned with the magnetic axis or with the line of sight, or it is slightly 
,. misaligned with respect to both the magnetic axis and the line of sight by 5 — 10° . 

' Conclusions. New observations, both from the ground and from the Hubble Space Telescope (HST), are important to characterise 
. y .' the optical/near- ultraviolet (UV) spectrum of RBS 1774, to better constrain the values of ro, d, and To and measure the source 
proper motion from which indirect constraints on the source distance can be inferred. 

Key words. Optical: stars; neutron stars: individual RBS 1774 



1. Introduction 

One of the most intriguing results of the all sky sur- 
vey performed by the Rontgen Satellite (ROSAT) has 
been the discovery of seven X-ray dim and radio-silent, 
middle-aged isolated neutron stars (a.k.a. XDINSs, see 
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Haberl 2007 and Turolla 2009, for the most recent re- 
views). XDINSs stand apart with respect to most X-ray 
emitting INSs due to their pure blackbody X-ray spectra 
{kT w 40 — 100 eV), produced from the cooling of the neu- 
tron star surface. Shallow X-ray pulsations (P = 3 — 12 
s), likely from large and hot polar caps, are observed 
for all but one of them. The measurement of the period 
derivative P in some XDINSs (see Kaplan & van Kerkwijk 
2009 and references therein) yielded, assuming magneto- 
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dipolar spin-down, ages of 1 — 2 Myrs and rotational 
energy losses £^ ^ 2 — 5 x 10'^° ergs s^^, too small to 
power detectable magnetospheric emission. Interestingly, 
XDINSs appear to be endowed with relatively high mag- 
netic fields, B « lO^'^-lO^* G, as inferred from spin-down 
measurements and the detection of broad spectral features 
which could be attributed to proton cyclotron scattering 
and/or electron transitions in H/He-like atoms. Besides 
the X-rays, XDINSs are only detected in the optical. In 
radio, stringent upper limits have been obtained at 820 
MHz (Kondatriev et al. 2009), while the detections of RX 
J1308.6-h2127 and IRXS J214303.7-H065419 at 111 MHz 
(e.g. Malofeev et al. 2007) is yet unconfirmed. XDINSs 
feature a > 5 optical excess with respect to the extrapola- 
tion of the X-ray spectrum (see, e.g. Mignani 2011) and, 
at least in the two best-studied cases (RX J0720.4— 3125 
and RX J1856.5— 3754), their optical fluxes closely follow 
a Rayleigh- Jeans (R-J) distribution, with a possible addi- 
tional power-law (PL) component for the former (Kaplan 
et al. 2003). Whether the XDINSs optical emission is pro- 
duced from regions of the star surface which are larger 
and cooler than those producing the X-rays (e.g. Pons et 
al. 2002) is debated. 

Recently, we observed IRXS J214303.7-f 065419 
(a.k.a. RBS 1774), the last of the ROSAT discovered 
XDINS (Zampieri et al. 2001) and possibly that with 
the highest inferred magnetic field (ss 10"'^^ G; Zane et 
al. 2Q05|3j with the Very Large Telescope ( VLT) and we 
identified its optical cotmterpart (B=27.4± 0.2; Zane et al. 
2008). Interestingly, the neutron star flux was found to be 
a factor « 35 above the optical extrapolation of the XMM- 
Newton X-ray spectrum (fcT ~ 104 eV; Nh 3.6 x 10^" 
cm~^; Zane et al. 2005), and a slightly larger optical excess 
was measured with the Large Binocular Telescope (LBT) 
by Schwope et al. (2009). So far, nothing is known about 
the optical/IR spectrum of RBS 1774, for which only the 
two B-band flux measurements of Zane et al. (2008) and 
Schwope et al. (2009) are available. Prior to our detection 
of the optical counterpart, the neutron star was observed 
both at optical and infrared (IR) wavelengths (Lo Curto 
et al. 2007; Rea et al. 2007; Posselt et al. 2009) but no 
counterpart was found. 

In order to constrain the optical spectrum of RBS 
1774, we performed new observations with the VXTin the 
R band to disentangle the slope of a flat PL continuum 
from that of a R-J distribution and, thus, to determine 
whether the optical emission is non-thermal or thermal, 
i.e. associated either with the magnetosphere or the neu- 
tron star surface. 

This paper is organised as follows: observations, data 
reduction and analysis are described in Sect. 2, while 
results are presented and discussed in Sect. 3 and. 
Conclusions follow. 



^ This value has been inferred from the observations of an ab- 
sorption feature in the X-ray spectrum of the source. However, 
the recent measurement of the pulsar spin-down (Kaplan & 
van Kerkwijk 2009) implies a magnetic field of ~ 2 x 10^^ G 



2. Observations and data analysis 

2.1. Observation description 

We obtained optical images of the RBS 1774 field with the 
VLT Antu telescope at the ESO Paranal observatory on 
September 19 and 20, 2009. Observations were performed 
in service mode with the FOcal Reducer/low dispersion 
Spectrograph {F0RS2), a multi-mode camera for imag- 
ing and long-slit/multi-object spectroscopy (Appenzeller 
et al. 1998), using the default RSpecial filter (A — 6550 
A; AA — 1650A). In order to achieve the highest sensitiv- 
ity at longer wavelengths, F0RS2 was equipped with its 
red-sensitive MIT detector, a mosaic of two 2kx4k CCDs 
optimised for wavelengths longer than 6000 A. In its stan- 
dard resolution mode, the detector has a pixel size of 0'.'25 
(2x2 binning) which corresponds to a projected field-of- 
view of 8f3 x 8f3 over the CCD mosaic. However, due to 
vignetting, the effective sky coverage of the two detectors 
is smaller than the projected detector field-of-view, and 
it is larger for the upper CCD chip. To include a larger 
number of reference stars for a precise image astrometry, 
as well as to increase the signal-to-noise ratio per pixel, we 
performed the observations in standard resolution mode. 
We positioned RBS 1774 in the upper CCD chip to exploit 
its larger effective sky coverage (7' x 4'). The instrument 
was operated in its standard low gain, fast read-out mode. 

To allow for cosmic ray removal, a sequence of 20 ex- 
posures of 775 s each were obtained in the two nights, for 
a total integration time of 9300 s and 6200 s for the first 
and second night, respectively. Exposures were taken in 
dark time and under photometric conditions, as recorded 
by the ESO ambient condition monitoJl and close to the 
zenith, with an airmass mostly below 1.3. The image qual- 
ity (IQ) was not constant during our observations. The 
seeing values measured by the Paranal differential image 
motion monitor (DIMM) are relative to the zenith and not 
to the pointing direction of the telescope and are subject 
to ground-layer turbulence (see also Martinez et al. 2010). 
Thus, they are not necessarily indicative of the actual IQ. 
For this reason, we computed it from the measured point 
spread function (PSF), derived by fitting the full width 
half maximum (FWHM) of a number of well-suited field 
stars using the Sextractor tool (Bertin & Arnouts 1996), 
following the recipe used by the F0RS2 data quality con- 
trol procedure^. The average IQ in the first night was 
■~ 0'.'6, while in the second night it was ^ 0'.'9. Bias, twi- 
light flat -fields frames, and images of standard star fields 
(Landolt 1992) were obtained as part of the F0RS2 sci- 
ence calibration plan. 

2.2. Data analysis 

Data were reduced through the ESO F0RS2 data re- 
duction pipelin43 for bias subtraction, and flat-field cor- 

^ http:/ /archive. eso.org/asm/ambient-server 
^ www. eso .org/observing/dfo/quality/FORS2/qc/qcl.html 
www.eso.org/observing/dfo/quality/FORS2/pipeline 
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rection. Photometric calibration was applied using the 
extinction-corrected night zero points computed by the 
F0RS2 pipeline and available through the instrument 
data quality control database. We converted the F0RS2 
zero points, computed in units of electrons/s, to units of 
ADU/s by applying the corresponding electrons-to-ADU 
conversion factors. We finally co-added the reduced sci- 
ence images using the IRAF task drizzle applying a Sct 
filter on the single pixel average to filter out residual hot 
and cold pixels. In order to achieve a better signal-to- 
noise, we decided to use only the images taken on the first 
night, which have the best IQ in our sample and provide 
the longest integration time (9300 s). Of course, we ver- 
ified that adding the images taken on the second night 
does not improve the signal-to-noise. 

To locate the RBS 1774 candidate counterpart in our 
F0RS2 image, we used as a reference its position mea- 
sured in our original FORSl observation. Since the FORSl 
detector was also a mosaic of two 2kx4k CCDs (but blue- 
optimised), with the same pixel size and field-of-view 
as F0RS2 at standard resolution, we decided to register 
the counterpart position on our co-added F0RS2 image 
through a relative astrometry procedure. Whenever pos- 
sible, this approach yields a position registration in the 
detector reference frame which is more accurate than that 
achievable through absolute astrometry, since it is not af- 
fected by the intrinsic absolute accuracy of the used refer- 
ence catalogue. Moreover, most of the reference stars from 
e.g. the GSC-2 (Lasker et al. 2008) are saturated in our 
images. For our relative astrometry, we used as a refer- 
ence the positions of several point-like objects detected in 
the field-of-view with Sextractor, after filtering out ob- 
jects that are either saturated, or too faint, or too close to 
the CCD edges. The relative frame registration was per- 
formed using standard tools available in ZR^Fand turned 
out to be accurate to better than 0.05 pixels (^ 0'.'0125) 
in both the x and y directions, which are aligned with 
right ascension and declination to better than 0.1°. To 
this, we added the uncertainty on the candidate counter- 
part centroid on the FORSl image which is ~ 0.1 pix- 
els (~ 0'.'025). Since the FORSl position refers to July 
11, 2007 (MJD=54292), we have to account for an addi- 
tional positional uncertainty due to the unknown neutron 
star proper motion to the date of our F0RS2 observations 
(MJD=55093). This corresponds to a yearly displacement 
of 0.2"tiioo/rfioo, where wioo and c?ioo are the unknown 
neutron star transverse velocity and distance^ in units of 
100 km s^^ and 100 pc, respectively. For, e.g. an average 
neutron star velocity of ^ 400 km s^^ (Hobbs et al. 2005) 
and a distance of 300 pc, this would translate into an ad- 
ditional positional uncertainty of ^ 0'.'6, which obviously 
dominates over the uncertainty on our relative astrometry. 



^ Posselt et al. (2009) quote a fiducial lower limit of 300 pc, 
estimated from models of the Nh distribution. 




Fig. 1. VLT/F0RS2 20" x 20" co-added R-band image 
(9300 s integration time) of the RBS 1774 field (North to 
the top, east to the left). The computed position of the 
candidate counterpart is marked by the cross. For a better 
visualisation, the cross arms (0"25) are equal to 20 times 
the la uncertainty of our relative astrometry. The circle 
(^ 0"6 radius) corresponds to the position uncertainty 
due to the unknown neutron star proper motion for, e.g. 
a transverse velocity of 400 km s~^ and a distance of 300 
pc. 

3. Results 

3.1. Photometry 

The computed position of the RBS 1774 candidate coun- 
terpart, overlaid on the co-added F0RS2 R-band image, 
is shown in Fig. 1. As it is apparent, no source can be 
visually identified at the expected neutron star position. 
We ran an automatic search with Sextractor but we did 
not detect any object above the 3a detection threshold. 
We applied an image smoothing using a Gaussian func- 
tion over 3x3 pixel cells but we could not detect any 
new feature which could be unambiguously associated to 
a real object and not to an enhanced background fiuctu- 
ation. Thus, we conclude that the optical counterpart of 
RBS 1774 is undetected in our images. We note that the 
0"6 positional uncertainty shown in Fig. 1 corresponds to 
the assumed values of the neutron star distance (300 pc) 
and transverse velocity (400 km s^^) and, thus, it is is 
only indicative. However, our conclusions are not affected 
by the actual position uncertainty of the neutron star since 
no new object is detected within a radius of 5" from the 
nominal candidate counterpart position with respect to 
those already detected (see Fig. 1 of Zanc et al. 2008). 
We determined the detection limit of our F0RS2 image 
from the number of counts corresponding to a 3a detec- 
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tion in an aperture of 2 pixel radius (O'.'S), derived from 
the standard deviation of the background sampled around 
the expected neutron star position. We then corrected the 
number of counts for the airmass, using the atmospheric 
extinction coefBcients measured for F0RS2 and available 
through the instrument data quality control database, and 
we applied the aperture correction. The latter was com- 
puted from the measured growth curve of a number of 
relatively bright but unsaturated stars in the field, with 
no adjacent star closer than 5", and whose photometry is 
not affected by strong background gradients. The final 3cr 
limit of 7? ^ 27 is about 1.3 magnitudes fainter than the 
r' upper limit reported by Rea et al. (2007). 

The non-detection of the Zane et al. (2008) candidate 
in our new VLT R-band images might cast some doubts 
as whether this object is indeed the optical counterpart 
of RBS 1774, although the low chance coincidence prob- 
ability with the Chandra position (2 x 10"'^; Zane et al. 
2008) already made the association quite robust. Given 
the flux of the candidate counterpart [B — 27.4 ±0.2), its 
non-detection in the R band corresponds to an observed 
(B — R) ^ 0.6. The observed upper limit on the colour, to- 
gether with the brightness of the candidate counterpart, 
can be used to rule out that this is a background ob- 
ject, positionally coincident with the neutron star. Since 
RBS 1774 coordinates point at a quite high galactic lat- 
itude {I = 62.6556°, b = -33.1392°), as shown by the 
many galaxies identified in the field, the candidate coun- 
terpart might be, in principle, an AGN. Indeed, most 
AGNs at z ^ 2 have ^ {B ~ R) ^ 3, which is com- 
patible with the colour of RBS 1774 candidate counter- 
part. The fact that no transient/persistent radio source 
was detected at its position in recent, deep 840 MHz ob- 
servations (Kondriatev et al. 2009) would only rule out 
a radio-loud AGN. However, an AGN not detected in X- 
ray^ would be either very much absorbed or at a much 
higher redshift, which would be hardly compatible with 
the constraint on its {B — R) and with its detection in the 
B band, respectively. On the other hand, a {B — R) ^ 0.6 
would be compatible with a main sequence star of spec- 
tral type earlier than F but this should be well outside 
the Galaxy to reproduce the observed B-band flux. Thus, 
both the counterpart colour and its brightness are only 
compatible with a neutron star. 

3.2. Multi-band spectrum 

We compared our i?-band magnitude upper limit with the 
best-fit model to the XMM- Newton spectrum of RBS 1774 
and with the available multi-band photometry (Rea et al. 
2007; Zane et al. 2008; Schwope et al. 2009). We com- 
puted the interstellar extinction correction from the hy- 
drogen column density Nh derived from the X-ray spec- 
tral fit. To this aim, we fitted the XMM-Newton spec- 
trum anew using updated calibrations and response files. 

^ Chandra observations of RBS 1774 rule out a blend with 
a background X-ray source. 



We re-extracted the EPIC-pn spectrum from the original 
May 2004 observation using version 10.0.0 of the XMM- 
Newton Science Analysis System (SAS), a circular region 
of 20" radius, and selecting only single-pixel events. Since 
combined fits with pn and MOS spectra are dominated 
by the higher statistical quality of the pn spectrum, we 
concentrated on the pn spectrum only. For comparison 
with the results obtained by Schwope et al. (2009) and 
Cropper et al. (2007) we fitted the spectrum with either a 
black body (BB) model, a BB plus an (additive) Gaussian 
line, and a BB plus a (multiplicative) absorption edge. In 
each case, interstellar absorption was included by using 
the XSPEC (version 12.6.0k) model phabs with abundances 
from Wilms et al. (2000). The best fit parameters are sum- 
marised in Table 1. As in previous analyses, inclusion of 
an absorption feature improves the fit over a simple BB 
model, although the current energy resolution does not al- 
low to discriminate a Gaussian line in absorption from an 
edge. The slight differences between the model parameters 
reported here and those given by Schwope et al. (2009) 
and Cropper et al. (2007) are largely due to the different 
revisions of the soft spectral response calibration of the 
EPIC-pn camera and possibly different elemental abun- 
dances used for the absorption model. The new spectral 
fit with an absorbed BB plus an (additive) Gaussian line 
yields a best-fit temperature of kT = 104.2 ± 2.1 eV, con- 
sistent with that obtained by Zane et al. (2005) by fitting 
the same XMM-Newton data but with older calibrations 
and response files. However, we note that our new spec- 
tral fit yields a best-fit Nh value of (2.6 ± 0.2) x 10^° 
cm^^ which is somewhat lower than that obtained by 
Zane et al. (2005), Nh = 3.6 x 10^° cm'^, and which was 
used both in Rea et al. (2007) and in Zane et al. (2008) 
to compute the interstellar extinction. This resulted in 
slightly larger extinction-corrected fluxes, where the effect 
is ^ 0.08 magnitudes, i.e. well below the uncertainties on 
their absolute photometry. Nonetheless, in the following 
we assume as a reference Nh = 2.6 x 10^" cm~^, obtained 
from our updated spectral fits to the XMM-Newton spec- 
trum. From this value, we derived an interstellar redden- 
ing E{B -V)^ 0.046 using the relation of Predehl & 
Schmitt (1995)0 and, from this, we computed the inter- 
stellar extinction in the different bands using the extinc- 
tion coefficients of Fitzpatrick (1999). We then corrected 
the available multi-band photometry accordingly. 

The new extinction-corrected multi-band fluxes from 
Rea et al. (2007), Zane et al. (2008), and Schwope et 
al. (2009) are shown in Fig. 2, together with the best-fit 
XMM-Newton spectrum and its extrapolation in the opti- 
cal domain. To these points, we added flux measurements 



^ This relation is affected by uncertainties for close objects, 
due to the problems of modelling the interstellar medium (ISM) 
at small distance from the Sun where microstructures weight 
more. We checked that, when using the relations of Bohlin et 
al. (1978) and of Paresce (1984) extinction corrections are con- 
sistent within 0.04 magnitudes, well below the pure statistical 
error on the source count rate. 
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Model 


Nh 


kT 


Eline 




X,.2(dof) 




(10^" cm- 2) 


(eV) 


(eV) 


(eV) 




(1) bbody 


2.7±0.2 


102.0±1.2 






1.51 (133) 


(2) bbody-|-gauss 


2.6±0.2 


104.2±2.1 


728±16 


31±31 


1.20 (130) 


(3) bbody*gauss 


2.6±0.2 


104.1±2.1 


733±17 


31±30 


1.20 (130) 



obtained in the near-uhraviolet (UV). The RBS 1774 field 
has been observed with the GALEX satellite (Martin et 
al. 2005) in both the NUV (A = 2771 A; AA ~ 530 A) and 
FUV passbands (A = 1528 A; AA ~ 220 A). We retrieved 
the fully reduced and calibrated imaging data from the 
GALEX archive[f| for inspection but we could not detect 
the source down to 3a upper limits of 2.76 ^Jy and 5.77 
/zJy in the NUV and FUV passbands, respectively. We 
corrected these fluxes for the interstellar extinction using 
as a reference the E{B — V) — 0.046 derived above and the 
extinction coefficients of Fitzpatrick (1999) for the NUV 
passband and of Seaton (1979) for FUV one. The field 
has been also observed with the XMM-Newton Optical 
Minotor (Mason et al. 2001) in the UWWl (A = 2675 
A; AA ~ 577 A), UWM2 (A = 2205 A; AA ~ 350 A), and 
UVW2 (A = 1894 A; AA - 330 A) filters but the source 
was not detected down to 3(T limits which are not deeper 
than the GALEX/NW one. At shorter wavelengths, the 
source was not detected by EUVE during its all-sky scan 
and was not targeted by pointed observations (Bowyer et 
al. 1996). As seen from Fig. 2, the GALEX points lie at 
least two orders of magnitude above the X-ray spectrum 
extrapolation. Thus, they can not be used to constrain the 
shape of the optical spectrum. 

3.3. The optical excess 

We note that the extrapolation of the new best-fit spec- 
trum in the optical domain lies a factor 1.35 above the 
extrapolation derived in Zane et al. (2008). Thus, as the 
result of the new fit to the XMM-Newton spectrum the 
optical excess of the VLT B-band flux is ^ 24.4 ± 3.7 (Icr 
confidence level). This is lower than the value of 35±7 (Icr) 
published in Zane et al. (2008) although still marginally 
consistent within la errors. As shown in Sect. 3.2, the dif- 
ference in the A'h obtained from our spectral fit and that 
of Zane et al. (2005) only corresponds to a difference of 
<j 0.08 magnitudes in the dereddened B-band fiux. Thus, 
it has only a negligible effect in the lower optical excess 
that we derive, which is mainly due to the difference in the 
extrapolation of the best-fit XMM-Newton spectrum, re- 
sulting from the use of new calibrations and response files. 
We also note that the value of the optical excess is insen- 
sitive to which of the two best-fit spectral models, (2) and 
(3) in Table 1, is actually assumed, since both the spectral 
parameters of the underlying blackbody and the inferred 



A^H are virtually identical. We note that Schwope et al. 
(2009), using LET data, obtained a slightly brighter B- 
band magnitude of 26.96±0.20 with respect to our value of 
B = 27.4±0.2 (Zane et al. 2008), which they interpreted as 
hint of a possible fiux variability from the source. Taking 
the magnitude difference at face value, this would imply 
a difference of « 50% in the optical excess. However, we 
note that the magnitude difference between the two mea- 
surements, 0.44 ±0.28 {la error), is not statistically com- 
pelling. Moreover, the two measurements have been taken 
with different telescope/detectors and with non-identical 
filters, which brings in an additional uncertainty due to 
the colour-term correction, and have been calibrated using 
different sets of standard stars, using default atmospheric 
extinction corrections. While Zane et al. (2008) applied 
the average B-band atmospheric extinction coefficient pro- 
vided by ESCd computed from observations performed at 
the Paranal Observatory over the April 2007-September 
2007 semester with the same telescope/instrument set-up 
used in their own observations, Schwope et al. (2009) ap- 
plied the B-band extinction coefficients measured at the 
LET site a few months apart, during the commissioning 
of the instrument. Since the measured atmospheric ex- 
tinction coefficients can display a night-to-night scatter 
up to ~ 0.1 magnitudes, this would introduce, on both 
measurements, an additional uncertainty of up to ^ 0.1 
magnitudes on the absolute photometry, for a target ob- 
served at the zenith, which is obviously larger if the target 
is observed at larger zenith angles. Accounting for all the 
above sources of uncertainties, the difference between the 
two B-band fiux measurements are significantly smoothed 
out, and we regard them as consistent with in the errors. 
Thus, in the following we use as a reference a value of the 
optical excess of ~ 24.4 ± 3.7, corresponding to our VLT 
B-band measurement. 

4. Discussion 

4.1. Non-thermal emission 

Our new R-band magnitude upper limit corresponds to an 
unabsorbed spectral flux < 1.33 x 10"^ keV (keV cm"^ 
keV^"'^), at a central energy of 1.77 eV. As it is seen from 
Fig. 2, the R-band upper limit is consistent with a set of 
PLs oc ly-"' with spectral indeces -2.07 < a < 0.07 
which pass ±3ct above/below the B-band fiux (shaded 
area). If the optical spectrum of RBS 1774 were non- 



www .galex.stsci.edu 



www.eso.org/ observing/ dfo / quality /FORS 1/qc/ qcl .html 
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log E (keV) 

Fig. 2. Optical/X-ray spectral energy distribution of RBS 1774. The solid line represents the updated best-fit model 
spectrum (fcT ^ 104 eV; Nh ^ 2.6 x 10^" cm~^) to the XMM-Newton data. The arrows represent the 5a optical/IR 
flux upper limits reported by Rea et al. (2007) and the 3a near/far- ultraviolet flux upper limits derived from archival 
G^ii?X observations (see text). The VXTB-band flux of the RBS 1774 candidate counterpart is marked by the cross, 
while our VLT R-band upper limit is marked by the cross and the arrow. The LBT B-band measurement (Schwope 
et al. 2009) is marked by the triangle. The shaded areas represent the region described by the two PLs with slopes 
a = —2.07 and 0.07 which pass ±3cr above/below the B-band flux and are consistent with our R-band upper limit (see 
text for details). Optical points have been corrected for the Nh obtained from the updated fit to the XMM-Newton 
spectrum (see text). The dashed line represents the blackbody at Tq = 11 eV (ro = 15 km, d =150 pc), normalized 
to our VLT B band point and compatible with the R-band upper limit. 



thermal, such a range of PLs would be only marginally 
compatible with the spectral indeces measured for other 
types of isolated neutron stars with non-thermal opti- 
cal emission, like the rotation-powered pulsars for which 
< a ^ 1 (see Mignani et al. 2007b, 2010a,b). However, 
since the actual R-band flux of RBS 1774 is obviously be- 
low the derived upper limit, the PL would most likely have 
a negative spectral index, which makes the spectrum rising 
and not declining towards the near-UV. Thus, although we 
can not rule out the presence of a PL component of non- 
thermal nature, we regard it as unlikely that the optical 
spectrum of RBS 1774 can be described by a single PL. 
Moreover, it is not clear what the origin of such a non- 
thermal PL could be. In Zane et al. (2008) we already 
deemed unlikely that non-thermal optical emission from 
RBS 1774 can be powered by the neutron star rotational 
energy loss which, if of the order of 10'^° erg as ob- 
served in the other XDINSs, would imply an anomalously 
large optical emission efficiency, ^ 1000 times larger than 
that of rotation-powered pulsars (Zharikov et al. 2006). 
Alternatively, on the basis of a detection in one band only 
and lacking a deep upper limit in the R band, we suggested 
that non-thermal optical emission could be powered by the 
presumably large neutron star magnetic field (« 10^^ G; 



Zane et al. 2005), as it has been proposed for the magne- 
tars, whose optical/IR emission efficiency is also a factor 
of w 1000 larger than that of the radio pulsars (Mignani 
et al. 2007a). However, in the few cases where multi-band 
photometry is available, magnetar optical spectra are con- 
sistent with PLs with spectral index a ^ 0, again very 
much different from our limit for RBS 1774. 

4.2. Thermal emission from the neutron star surface 

The R-band upper limit is consistent with a R-J spectrum 
which passes through the B-band flux. A linear function 
connecting the VLT B-band point and the R-band upper 
limit has a slope ^ —1.26, and the R-band upper limit lies 
only ^ 25% above the extrapolation of the R-J spectrum 
that passes through the central value of the measured B- 
band flux. This is well compatible with both the B and R- 
band fluxes being on the same R-J spectrum, conceivably 
emitted by a region of the star surface which is colder than 
that responsible for the X-ray emission. 

As in Zane et al. (2008), we consider a picture in which 
the optical emission originates from a colder fraction of the 
neutron star surface, which emits a blackbody spectrum 
at a temperature Tq (Braje & Romani 2002; Pons et al. 
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X (deg) X (deg) 

Fig. 3. Left: the contour of constant PF=0.04 (solid line) for the case ro = 15 km, d — 150 pc, Tq = 11 eV, rx = 1 
km, and M = 1.4 Af0. The dashed line shows the analytical result for point-like caps (see text for details). Pulse 
profiles are of class I according to the classification of Beloborodov (2002), i.e. they do not exhibit flat intervals, only 
when the contour is below the straight line. Right: same for M — 1.8 Mq. 



2002; Kaplan et al. 2003; Triimper et al. 2004). In the 
R-J tail, the ratio between the optical and X-ray fluxes 
scales as ~ r^To/r^Tx = f, where ro is the size of the 
optically emitting region (which, of course, cannot exceed 
the value of the neutron star radius), while rx and Tx 
are the radius and temperature of the X-ray emitting re- 
gion, as inferred from X-ray spectroscopy, Tx = 104 eV, 
rx = 2((i/300pc) km. Since the RBS 1774 distance is un- 
known, we iterated our computations over several values 
of d and ro and derived the value of To that corresponds 
to a given value of the optical excess, / (as computed 
from the B-band flux) . Since no contribution from such a 
cold component is observed in the 0.1-1 keV XMM-Newton 
spectruno, it must be i? <C 1, where 

^ /roVfToV Io.% t'/bMt) - 1]^^ 

Furthermore, the extrapolation of the cold blackbody 
component in the R band must be compatible with our 
new VXT upper limit. We have performed the calculation 
for a large set of parameters, varying the source distance 
between 100 and 500 pc, the B-band excess / within the 
3(7 limits (13.3 < / < 35.5), and considering ro = 10, 
15, and 20 knS We found that the data are quite con- 
straining: a radius ro = 15 km is compatible with the 
considered range of / only for d < 200 pc (which implies 

A soft thermal component to the X-ray spectrum could be 
compatible with the XMM-Newton RGS data which show an 
emission excess at low energies (Schwope et al. 2009) 

Neutron star radii of ^ 20 km are unlikely for most equa- 
tions of state (e.g., Lattimer & Prakash 2004). 



fx < 1-3 km) and To < 20 eV. Assuming a smaller ra- 
dius (ro = 10 km), the observed excess is only compatible 
with even smaller distances (^ 100 pc, corresponding to 
rx = 0.6 km) and Tq < 10 eV. On the other hand, assum- 
ing a larger radius (ro = 20 km) the observed excess is 
compatible with distances up to 300 pc (corresponding to 
fx ^ 2 km) and Tq < 15 eV. As an example, we plotted 
in Fig. 2 a blackbody with To = 11 eV, corresponding to 
ro — 15 km and d = 150 pc. 

Next, we verifled if the inferred size of the hotter caps 
is compatible with the observed pulsed fraction (PF) 

PF ~ — Fmin) ^ q ^2) 

where Fmax {Fmin) are the maximum (minimum) value 
of the flux along the pulse The PF was computed nu- 
merically by means of an IDL script for different allowed 
combinations of ro, rx and To (that is for which the 
condition R 1 was met), and assuming in all cases 
Tx — 104 eV. We also used two values of the star mass, 
Al — 1.4, 1.8 Mq to check how sensitive our results are 
to general relativistic effects (which depend on the com- 
pactness M/R; e.g Beloborodov 2002). Since it is always 
rx ^ O.lro, we took the latter to coincide with the star 
radius. Results are shown in Fig. 3, where the curve of 
constant PF (0.04) is plotted as a function of the two 
geometrical angles, f , the angle between the neutron star 
spin and magnetic axes, and x, the angle between the spin 
axis and the line-of-sight (LOS), for the case reported in 
Fig. 2. Overplotted to the numerical contour is the ana- 
lytical curve obtained following the method presented by 
Beloborodov (2002), and which is strictly valid only for 
point-like caps. The present analyses confirms previous 
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findings (Zane et al. 2008; Schwope et al. 2009). One pos- 
sibility is that either the spin axis is nearly aligned with 
the magnetic axis, hence the star is a nearly aligned ro- 
tator, or the spin axis is nearly aligned with the LOS. 
The other possibility is that the spin axis is slightly mis- 
aligned with respect to both the magnetic axis and the 
LOS by 5 — 10°. We note that although the numerical and 
analytical curves are qualitatively similar there are quan- 
titative differences. The allowed range in the two angles 
is somehow smaller in the complete calculation, as it is 
expected when the finite (albeit small) size of the caps 
is accounted for. This also may explain the somehow dif- 
ferent conclusions reached by Schwope et al. (2009), who 
used a treatment valid only for point-like caps. The small 
opening of the X-ray emitting regions does not constrain 
the geometry too much, and several combinations of the 
angles ^ and x &re possible. 

5. Summary and conclusions 

We performed deep optical observations of RBS 1774 in 
the R band with F0RS2 at the VLT to obtain the first 
characterisation of its optical spectrum. We did not de- 
tect RBS 1774 candidate counterpart at its expected po- 
sition down to a 3cr limiting magnitude of R ~ 27. The 
constraint on the colour of the candidate counterpart, 
{B — R) ^ 0.6, rules out a foreground object, position- 
ally coincident with the X-ray source. We re-analysed the 
XMM-Newton data of RBS 1774 using new calibrations 
and response files (see Sect. 3.2). We found that the op- 
tical excess of the VLT B-band flux (Zane et al. 2008) 
with respect to the best-fit X-ray spectrum extrapolation 
is now 24.4 ± 11.1 (3cr confidence level). This value is still 
incompatible with rotation-powered non-thermal emission 
from RBS 1774, unless its optical emission efficiency is 
~ 1000 times larger than that of rotation-powered pul- 
sars (Zharikov et al. 2006). Moreover, our R-band upper 
limit would most likely imply a PL spectral index a < 0, 
while rotation-powered pulsars have usually < a < 1 
(Mignani et al. 2007b; 2010a,b). On the other hand, ex- 
plaining the optical excess in terms of pure thermal emis- 
sion from the neutron star surface would require both a 
small distance and rather stringent limits on the inclina- 
tion of the LOS and the magnetic axis with respect to the 
neutron star spin axis. For instance, a radius rp < 15 km 
implies a neutron star distance d < 200 pc, which im- 
plies To < 20 eV and rx < 1.3 km. Such a small X- 
ray emitting area is compatible with the ~ 0.04 X-ray 
pulsed fraction if either the star spin axis is nearly aligned 
with the magnetic axis or with the LOS, or it is slightly 
misaligned with respect to both the magnetic axis and 
the LOS by 5 — 10°. Observations in the near-UV would 
obviously be crucial to constrain the slope of the RBS 
1774 optical spectrum and to better constrain Tq and rp- 
With the GALEX fluxes falling well above all possible R- 
J spectra compatible with the VLT flux measurements, 
only the Hubble Space Telescope (HST) can provide the 
required near-UV sensitivity. Future optical observations. 



both from the ground and from space, will also be im- 
portant to measure the source proper motion, to inde- 
pendently confirm the optical identification and to infer 
indirect constraints on the distance. 

Acknowledgements. RPM thanks the ESO User Support 
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